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I. INTRODUCTION 



Until recently, few al 
hydraulic or pneumatic act 
applications. However, a 
earth, permanent magnet mat 
ductor transistor technolo 
electromechanical actuator 
Utilizing rare earth (s 
magnetic material within 
provides a prime mover 
offering superior performan 
ydraulic systems. 

Elimination of the brus 
the dc motor provides numer 
motor speed along with a re 
a given horsepower, 2) th 
rotors instead of a rotatin 
the elimination of the 
design, implementation and 
since there are no brushes 
allowing motor operation in 
improved thermal characte 
core) , which arise primar 
of the machine, are easil 
housing. 

Use of the brushless d 
cerain disadvantages , chie 
the uncertain availability 
for rotor constr uction. 
conducted in which the perf 
ferrite and a samarium 



ternatives existed for the use of 
uators in primary flight control 
dvances made in the field of rare 
erials and in high power semicon- 
gy has made possible the use of 
s as a practical alternative, 
pecifically, samarium cobalt) 
a brushless dc motor design 
in flight control applications 
ce characteristics over electroh- 

h type commutation scheme within 
ous advantages: 1) higher rated 

duction in weight and volume for 
e ability to use permanent magnet 
g armature winding, combined with 
brush assembly translates into 
maintenance simplifications, 3) 

, no arcing will occur, hence 
hazardous environments, and 4) 
ristics, as losses (ohmic and 
ily within the stationary portion 
y dissipated through the stator 

c machines does not come without 
f among them being the cost and 
of the samarium cobalt material 
As an aside, a recent study was 
ormance characteristics of both a 
cobalt type dc motor w«=re 
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investigated [Ref. 1]. The research conducted demonstrated 
the superiority of the samarium cobalt design, and hence its 
desirability for present and future applications. An 
additional disadvantage accrues from the fact that motor 
commutation must be accomplished electronically, resulting 
in an increase in the complexity of the design of the motor 
controller, with a commensurate increase in the cost of its 
implementation . 

This study will concentrate on the power control of 
brushless dc motors utilizing a pulsed power approach. Power 
pulse control of dc motors, better Jcnown as pulsewidth modu- 
lation (EWM), utilizes as an input to the motor, voltage or 
current pulses with the pulse duration being controlled. 
Pulsewidth modulation offers considerable advantages in the 
control of dc motors, which will be outlined in the 
following chapter. The research was conducted utilizing a 
computer model of a 3-phase, 4 pole brushless dc motor, 
which was developed by Thomas in a related study [Ref. 2]. 

It is to be noted that this research project represents 
only one part of an ongoing effort to accurately simulate 
the use of a brushless dc motor as an electromagnetic actu- 
ator for use in advanced missile control systems. While 
electromagnetic actuators have previously seen use in 
missile projects such as HARM and Condor, these actuators 
have been too large and have had too slow a response for 
high torque applications as are found in the AMfiAAM 
(Advanced Medium Range Air to Air Missile) project [Ref. 3]. 
The overall aim of the entire research program is to help 
exploit the technological opportunity which exists within 
rare earth permanent magnet dc machines in the role of elec- 
tromechanical actuators for use in guided missile control 
appl ications. 
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II- POLSEWIDTH MODULATION 

Pulsewidth modulated switching amplifiers offer consid- 
erable advantages in the control of dc motors. Before 
describing the techniques used in simulating the pulsewidth 
modulated control of the computer model mentioned previ- 
ously, it will be important first to outline the basic prin- 
ciples and characteristics of pulsewidth modulation. 

A. PULSEWIDTH MODULATION PRINCIPLES 

The pulsewidth modulation scheme utilizes transistors in 
the switching mode, whereby the transistors are switched 
into and out of saturation. This switching action results 
in the minimization cf power losses in the transistors, with 
a savings realized in reduced heat sinking requirements and 
in the usage of less expensive power transistors. Since the 
power transistors are switched on and off at a frequency 
beyond the system bandwidth, the motor will filter the high 
frequency components of the modulated signal and respond 
only to the signal’s low frequency components. 

Closing a feedback loop around the pulsewidth modulated 
amplifier results in the amplifier behaving as either a 
current or a voltage source. The feedback loop allows the 
motor to be connected with the amplifier through additional 
series inductance, resulting in a smoothing of any current 
ripple. The feedback loop also allows for easy current 
limiting, merely by limiting the output of the feedback 
summing amplifier. Finally, a feedback loop provides output 
short circuit protection, as the output current is deter- 
mined by the input voltage without regards to the output 
impedance- 
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B. MODULATION TECHNIQUES 



There are two basic methods for obtaining a pulse width 
modulated signal. 

1 • The Dit h er Hethod 



The first technique requires that the input signal 
(Xo ) f be added to a high frequency sawtooth signal (also 
known as a dither signal). After summing, the resultant 
signal (Y1) , is fed into a relay element. The relay element 
converts the summed signal into a two level output (Y2) , 
which then switches from +V to -V whenever Y1 experiences a 
zero crossing, as shown in Figure 2.1. 

The duty cycle (a) of the output signal is related 
to the input signal and magnitude (E) of the sawtooth signal 
by: 



a = (E + X)/2E 



(eqn 2.1) 



Utilizing this technigue has as an advantage the fact that 
one may control the frequency of the supplied sawtooth 
signal without any changes to the motor controller 

circuitry. 

The Limit Cycle Method 

The second technique for producing a pulsewidth 
modulated signal is by closing a feedback loop around a two 
level switch. The feedback signal causes the system to 
exhibit high frequency limit cycle behavior. Taking a 
velocity control system as an example, if the velocity error 
signal (referance signal minus the actual velocity) were 
negative, indicating that the velocity was too high, the 
output signal from the switch would tnen serve to drive the 
velocity lower. As the velocity of the system dropped below 
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time 






Figure 2.1 Creation of a PWM Signal 
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the reference level, the 
positive, causing the oupu 
and consequently force the 
In this manner the switch o 
tions, describing a pulsetr 
by the voltage levels at w 
the dynamics in the feedba 
velocity waveform would ne 
the output of the switch, 
would become constant and p 
signal at the output of 
specific system dynamics, 
operating under load, the 
the lead on the motor, 
cycle velocity control sche 



error signal would then become 
t of the switch. to ohange states 
velocity to once again increase, 
utput would thus exhibit oscilla- 
ain whose frequency is determined 
hich tae switch operates and by 
ck path. The frequency of the 
cessarily be the same as that of 
and under steady state conditions 
eriodic. The pulsewidth of the 
the switch is dependent upon 
For tae case of a d.c. motor 
pulsewidth is directly related to 
The specific details of a limit 
me is examined in Chapter 3. 



C. ANALYSIS OF THE PULSEWIDTH MODULATED SIGNAL 



A diagram of an ideal 
The "freewheeling" diode ( 
during the pulse off period 
to circulate. Figure 2.3 
current and voltage relat 
scenario. 

Because the supply volt 
cies typically on the orde 
important to study the pow 
the power being pulsed on a 
the evaluation of power 1 
armature resistance is seen 



dc motor is shown in figure 2.2. 
FWD) serves to bypass the motor 
, allowing the armature current 
shows the typical steady state 
ionship in a pulsewidth control 
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Figure 


2. 2 


Bas ic 


DC Motor Schematic 


We 


ma y 


now 


define 


the 


curce nt 


form factor (k) as the ratio 


of 


the 


RMS 


current 


to 


the aver 


age current: 



“ = ^rms/^ave 2.3) 

Developed Eotor torgue is directly proportional to motor 
current (equation 2.4), allowing motor losses (under ?KM 
conditions) to be described as in equation 2.5, 

Tg = (©9^ 2.4) 

?! = R*k2*I2gyg (egn 2.5) 

Substituting equation 2.4 into equation 2.5 demonstrates 
that motor losses are dependent on the current form factor 
and the output motor torque, as seen below: 

= (E/K_^2) *k 2*T^2 (eqn 2.6) 
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1 




Figure 2.3 PHM Steady State Behavior 



The IDO t or forni factor has a large influence on niotor 
heating and hence, power losses. Since performance in speed 
control systems is often limited by power dissipation 
constraints, it is important to determine the form factor 
for a given PWM scheme. The relationship between motor 
armature losses and form factor is shown in Figure 2.4 
[Bef. 4]. 

In order to determine the form factor for a given 
system, one must first be able to determine a system's 
average and RHS currents. The differential equations 
describing the motor action for the basic dc motor system 
(as seen in Figure 2.2) are as follows: 



Pulse on: Ldi/dt = V - Rl - k^*w 

D 



(eqn 2.7) 
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Figure 2.4 Additional Armature Losses vs. Form Factor 

Pulse off: Ldi/dt = HI - Kb*w (egn 2.8) 

where Kt is the counter e mf constant and w is the motor 
speed. Solving the differential equation for motor current 
(and referring to Figure 2. 3) yields: 

Pulse on: I = I - (I - I )exp(-t/t ) (egn 2.9) 

n n oo ^ n ^ 

Pulse off: I = (Iq ♦ (egn 2.10) 
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where tn is the system electrical time constant, 



= (V - K^*w)/E (eqn 2.11) 

and 

= K^*w/E (eqn 2. 12) 

The average and EMS currents are then: 

I = (al - bl ) /F (egn 2. 13) 

ave ' n m ^ j 



12 = (Vp) ((al 2 + bl 2)- - t (I ) (I )) (egn 2,14) 

rms n m n oa t 

where 



a = t *ln({I - I )/(I - I )) (egn 2,15) 

n n oo n o 



b = t *ln ((I + I )/{I + I ) ) 

n o m oo m 



(eqn 2. 16) 



od 



= - IoO> 



(eqn 2. 17) 



I 



t 




+ I ) 

m 



(eqn 2. 18) 



and 1/P is the switching frequency. 



As a means of demonstra tin j the viability of PWM control 
of dc motors, simulations were conducted using fixed duty 
cycle power pulses to determine the relationship between 
form factor and motor load torque. A computer program was 
written to analyze the output motor current waveshape for 
the average and rms currents utilizing relationships 
detailed in this section. The program used to compute these 
currents may be found in Appendix C. Figures 2.5 and 2.6 
show that the form factor rapidly approaches unity as the 
load on the motor is increased, indicating that the motor is 
experiencing only slight additional losses (in terms of 
percentages) due to the power pulsing effect as compared 
with a constant voltage supply arrangement. It is also 

pointed out that the motor runs more efficiently (lower form 
factor) at higher frequencies for a given load torque. This 
is due to the fact that as the motor is pulsed more 
frequently, the motor speed will not drop off as fast and 
hence, the energy required to move the mass of the rotor 
back to its steady state speed will not be as great. 
However, switching losses in the transistors will usually 
limit the switching frequency to less than ten KHz. 



20 



FORM FACTOR 




Figure 2.5 Form Factor vs. Load Torque (1 KHz) 
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Figure 2.6 For» Factor vs. Load Torque (5 KHz) 
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III. PULSE WIDT H HODOLATED SPEED CONTROL 

The ability of a dc motor to maintain a given speed when 
a load torque is applied is generally referred to as speed 
regulation. Although a dc motor by itself is an open loop 
system/ the prescence of the back electromotive force (bemf) 
signal serves to close a natural/ "built-in"/ feedback loop, 
as shown in Figure 3.1 However, because the dc motor is 
intrinsically an open loop system with relatively constant 
power input, as the load torque increases, the speed will 
decrease, and hence, no speed regulation may be achieved. 
In order to maintain a constant speed, the input power must 
increase with the applied load. 




Figure 3. 1 Block Diagram of a Basic DC Motor 



Figure 
pulse width 



3.2 illustrates the 
modulated amplifiers. 



basic principle behind 
What is important to note 
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Figure 3.2 Pulse Width as a Function of Load Torque 
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A. COEfiENT FEEDBACK 

Prior to developing 
studies vere conducted of m 
pulsewidth power pulsed in 
Motor speed was found, as 
voltage, to vary linearly 
Figure 3.3. As an aside, 
the load torques ranging f 
as this range represented t 
the motor modelled in the 
of average current vs. loa 
simulations are shown in Fi 
• Studying the curves fou 
the conclusion that a schem 
could be developed with th 
proportional to the average 
of the motor pulsewidth was 

PW = DCF + * K 

O.V 0 

The motor parameters found 
follows : 

PW = the input pulse wid 
= average motor curr 
DCF = a dc term which es 
K = the current feedba 

It was decided that the 
50% duty cycle pulses to m 
occur at smaller duty cycle 
a frequency set at 5 KHz, 
speed at approximately 1375 
sients could be expected wh 



any specific control schemes, 
otor behavior relative to a fixed 
put, at a frequency of 5 KHz. 
in the case of constant supply 
with load torque , as shown in 
all studies were conducted with 
rom zero to eighty ounce-inches, 
he linear range of operation for 
Thomas study [Ref. 2]. The plot 
d torque for the fixed pulsewidth 
gure 3.4. 

nd in Figures 3.3 and 3.4 led to 
e for speed control of the motor 
e pulsewidth being made directly 
motor current. The basic form 
decided to be as follows: 

(eqn 3.1) 

in equation 3. 1 may be defined as 

th (duty cycle) 
ent 

tablishes no load speed 
ck constant 

motor would run at a minimum of 
inimize power losses which would 
s in light load conditions. With 
this necessarily fixed no load 
rpm. Since large current tran- 
en the motor was started or when 
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Figure 3.3 Motor Speed vs. load Torque (fixed pulse width) 
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Figure 3,4 Current vs. Load Torque (fixed pulse width) 
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the load condition was changed, a limiter was written into 
the simulation program to set the maximum pulsewidth at just 
less than the 100% duty cycle point. It was not set at 100% 
due to difficulties encountered with the simulation 

language. Other additions to the basic computer model 
included a limiter to prevent pulsewidths with less than 50% 
duty cycles as well as a current limiter to prevent negative 
currents. The latter was added to simulate the effect of 
the addition of the "freewheeling" diode described in the 
preceding chapter. 

To set the no load speed at 1375, it was necessary to 
obtain parameters that would establish PW equal to .5 at 
zero load torque. The value of K was determined from 
Figures 3.3 and 3.4 by noting that the 75% duty cycle condi- 
tion, at approximately 1350 rpm, occurred at a load torque 
of 40.0 oz-ins, which also corresponded to an average motor 
current of 2.50 amperes. Setting PR to .75 and DCF to .48 
in eqn. 3.1 led to a K value of .107. A value of DCF of .48 
results in an approximately 50% pulsewidth modulated signal, 
as no load motor current is approximately .002 amps. It was 
felt that higher speeds could then be achieved by increasing 
DCF, as the basic relationship between pulsewidth and speed 
appeared linear for any given load torque. 

Extensive simulations were conducted utilizing the feed- 
back control relationship as shown in equation 3.2. Graphic 
results for these simulation trials are shown in Figure 3.5. 

PW = .48 + (I^^^*.107) (eqn 3.2) 

It is clearly evident from Figure 3.5 that the control 
scheme utilized performed u nsatisfactorilly for its task of 
maintaining constant speed throughout the given range of 
load torques. While certain variations might have been 
expected in the output speed, the results demonstrated 
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nonlinearities in motor performance wnich were clearly 

unsuitable for its given application. One of the major 

problems encountered with this control scheme stemmed from 
the fact that the average current was used as a feedback 
parameter, rather than the actual motor current (the ripple 
present in the motor current was deemed to be too high to be 

used in a velocity control scheme based on current feed- 

back) - The average motor current proved to be unsatisfac- 
tory for the given task for two reasons: 1) average motor 

current requires time to approach the system’s actual 
average current value due to the changes to the current 
incurred by transients such as are caused by changes to the 
system’s steady state behavior and thus adds significantly 
to the motor’s settling time to variances in load or 
commanded speed, and 2) apparent non linearities which would 
appear if pulse duty cycle were plotted as a function of 
load torgue, which runs contrary to the initial assumptions 
upon which the control scheme was devised. 
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attempted 
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position error signal 
lation of the complete 
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Figure 3.5 Motor Speed vs. Load Torque (current feedback) 
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set point, the power is turned off. The controlled velocity 
waveshape is diagrammed in Figure 3.7. 




Figure 3.6 Schematic Diagram Of Pulsewidth Modulator 



Operation of the system 
straightforward. Referrin 
speed is below (VCOM + VTOL 
allowing motor current to 
(VCOM + VTOL), QM is swit 
large voltage across the 
rapid rate of change in the 
voltage turns on diode DM, 
decaying motor current. W 
- VTOL) , QM is switched ba 
diode DB are utilized whe 
regenerative braking mode o 
dc motor is used as a gener 



in the limit mode is relatively 
g to Figure 3.7, if the motor 
) , transistor QM is switched on, 
flow. When the speed reaches 
ched off, which then induces a 
inductor terminals, owing to a 
inductor current. This induced 
which provides a path for the 
hen motor speed decays past (VCOM 
ck on again. Transistor QB and 
n the motor is operated in the 
f operation. In this mode, the 
ator as the inertia of the rotor 
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Figure 3-7 Motor 

is converted into electrical 
dc voltage supply while con 
rotating masses. Again, ho 
function of the design of 
actuator, and will not be f 
control of transistors QB 
utili 2 ing voltage comparators, 
proportional to the commanded 
tolerance, while the other in 
to the actual motor speed, 
realize now is that the pulse 
are both variable, and will be 
of the system. 
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Velocity Waveshape 

energy to charge the system’s 
trolling the braking of the 
w this mode is utilized is a 
the motor's electromechanical 
urther commented upon. The 
and QM can be accomplished 
where one input is a voltage 
speed plus the allowed speed 
put is a voltage proportional 
What becomes important to 
width and the pulse frequency 
dependent on certain dynamics 



A trial simulation was conducted using the limit control 
scheme. The speed tolerance was set at five rpm, and the 
commanded speed was set to 1400 rpm. Table I contains the 
results of the simulation. This data is also represented 

graphically in Figure 3.8. Loading of the motor was accom- 

plished using a terminated ramp signal; the terminal value 
of the ramp is the desired motor loading. The speed accu- 
racy of the motor is defined as the difference between the 
maximum and minimum motor speed divided by the commanded 
speed . 



TABLE I 

Motor Characteristics 



Load Torque 
( 02 -in) 



maximum 



Speed (rpm) 

minimum 



0.0 1416.0 1385.3 

16.0 1418.2 1375.9 

32.0 1417.9 1369.7 

48.0 1416.1 1366.2 

64.0 1414.3 1362.1 

80.0 1409.6 1355.8 



% accuracy 



2.19 

3.00 

3.44 

3.57 

3.71 

3.79 



The results of the initial simulations indicates that 
positive control of the commanded motor speed may be 
accomplished utilising the limit cycle method. All further 
studies are therefore based on a computer model whose speed 
is controlled in this manner. 

Studies of system charac teristics (pulsewidth, pulse 
frequency, speed regulation) are presented in detail in the 
next chapter in order to better define the operational 
enevelope of the modelled dc motor using the velocity limit 
control scheme. Addittion ally, the effects of the addition 
of series inductance are also investigated to determine its 
affects in regards to ripple reduction. 
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5 SPEED ACCURACY 




Figure 3.8 Hotor Speed Accuracy vs. Load Torque 
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IV 



PERFORMANCE OPTIMIZATION 



The simulations described at the end of the previous 
chapter were not necessarily indicative of the optimum 
performance potential of the modelled dc motor. It appears 
evident that there must exist means by which the current 
ripple and the speed regulation may be respectively reduced 
and improved. One way to achieve improved speed accuracy 
(which implies that the ripple present in the motor velocity 
is reduced) is to decrease the speed tolerances established 
for motor operation. The effects of the variation of the 
speed tclearnce settings will be examined shortly. However, 
to reduce the current ripple, which in turn translates to 
reduced power losses in the motor, it has been suggested 
that one must add inductance in series with the motor arma- 
ture [Ref. 5]. The effect of the addition of series induc- 
tance on motor performance is studied in the following 
section. 



A. ADDITION OF SERIES INDUCTANCE 



A dimensionless current 
curent ripple multiplied by 
and divided by the supply 
itself is defined as the diff 
the time the motor is pulsed o 
is pulsed off, or simply a 
minimum and maximum currents, 

I = i(aT) - i(0) 



be defined as the 
armature resistance 
The current ripple 
erence between the current at 
n and the time when the motor 
s the difference between the 
as follows: 

(eqn 4.1) 



ripple may 
the motor 
voltage. 
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For th€ case where only one power supply is used (unidirec- 
tional drive), dimensionless current ripple 
defined as follows: 

= (1+exp(-t^) -exp(-a*t^)-exp ((1-a) tj^) )/g (eqn 4.2) 

where : 

t^ = L/R (eqn 4,3) 

and 

g = 1 - exp(-tj^) (eqn 4.4) 

Dimensionless current ripple may be plotted versus the 
pulsewidth, or duty cycle, with the ratio of the period of 
the PWM signal to the electrical time constant, tau, to form 
a family of curves, as in Figure 4.1. 

As can be seen from Figure 4.1, the magnitude of the 
current ripple depends to a great extent upon the ratio of 
the pulse period to the motor electrical time constant. 1-7e 
should therefore expect that that reduction of the motor 
current ripple can be accomplished through the addition of 
series inductance, which reduces the magintude of tau. 

The effect of adding series inductance to the motor may 
be seen by writing the differential equations for the 
simplified circuit diagram of a dc motor as shown in Figure 
4.2. Applying Kirchoff's Law and summing voltages around 
the Icop results in a current -voltage relationship as shown 
in equation 4.5: 

V = (L + L1)di/dt + Ri + K^*w (egn 4.5) 

Taking the Laplace transform of equation 4.5 yields: 
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DIMENSIONLESS CURRENT RIPPLE 
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Figure 4.1 Current Ripple vrs. Pulse Duty Cycle 
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V (s) = s (L +L1)I(sJ + RI (s) + Kj^-W(s) 



(egn 4.6) 



Solving equation 4.6 for motor current yields: 

I(s) = (l/s(L + LI) + R)).(V(s) - K^VHs)) 

From equation 4.7 the electrical time constant is: 

t^ = (L + LI) /R (eqn 4.8) 




Figure 4.2 Basic DC Motor Circuit Diagram 

It is tc be noted that the previous derivation ignores 
the affect of the resistance which will accompany the 
additional series inductance. However, it will be assumed 
that a value for the electrical time constant may be set as 
in equation 4.8 through careful selection of the type and 
quantity of series inductance added to the motor. 
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To test the effects of the addition of series induc- 
tance, the computer model was modified to reflect the addi- 
tion of inductance equal in magnitude to the inductance 
present in the motor's windings (ignoring the change in the 
total resistance for the reasons mentioned above). 
Simulations were then conducted throughout the load torque 
range of from zerc to eighty ounce-inches with motor 
velocity set at at 1400 rpm (approximately half of the 
motor's no load speed) and speed tolerance set to ±1 rpm. 
Table II summarizes the results of these simulation trials- 
The data for simulations made without the additional induc- 
tance is shown in Table III and is included for sake of 
comparison with the data in Table II. 



TABLE II 

Inductance Effects on Motor Operation 



Load Torque 
(oz-ins) 


la ve 
(amps) 


Irms 
(amp s) 


Form Factor 


Current 

Ripple 


05. 0 


0. 257 


0.454 


1.76 5 


. 1100 


16.0 


0.9 95 


1. 131 


1. 138 


.0874 


32. 0 


1 . 967 


2.032 


1. 033 


.0948 


48.0 


3.0 23 


3. 064 


1.012 


.0562 


64. 0 


4. 0 09 


4.028 


1. 005 


.0787 


80. 0 


5.011 


5.031 


1. 004 


. 1295 



TABLE III 

Motor Performance (no series inductance) 



Load Torque 
(oz-ins) 

05.0 
16. 0 

32.0 

48.0 
64. 0 
80. 0 



xa ve 
(amps) 


Irms 
(amp s) 


Form Factor 


Curre 

Rippl 


0. 179 


0.347 


1.939 


.326 


1.019 


1.10 1 


1.081 


.259 


1. 983 


2.026 


1.022 


.279 


3. 013 


3.370 


1.008 


. 261 


4. 028 


4.04 0 


1.003 


. 226 


5. 017 


5.026 


1.002 


.285 
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DIMENSIONLESS CURRENT RIPPLE 
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DIMENSIONLESS CURRENT RIPPLE 




Figure 4.3 Inductance Effects on Current Ripple 
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Examination of the data in Tables II and III demon- 
strates clearly the ripple reduction obtained through the 
addition of the series inductance. Motor form factor was 
calculated to show that power losses (at low load) in the 
transistors are also reduced, and tnus represents an addi- 
tional benefit gained from the added inductance. However, 
while series inductance will reduce the current ripple 
without affecting steady state behavior, it will have an 
affect on the transient behavior of the system, as will now 
be shown. 

Assuming now that the system is operating in steady 
state, a step input command (such as a change in the motor’s 
commanded velocity) will force the pulsewidth modulated 
signal to the full on condition. The response of the system 
will now be limited by the motor’s electrical time constant. 
In the case where series inductance has been added to the 

system, the response of the motor to a step input will be 

slower than if the additional series inductance were not 
present . 

To test this effect, simulation trials were conducted 
where the motor was allowed to achieve a steady state speed 
of 1000 rpm and then was subjected a step input command to 
increase motor speed to 1400 rpm. Torgue load was 32 oz-ins 
and the speed tolerance setting was ± 1. The response time, 

or the time required for the motor to settle at the new 

commanded speed, was measured for trials in the motor’s 
standard configuration and for the case where series induc- 
tance was added. For the first case, where there was no 
additional inductance, the response time was measured at 
approximately 2.5 milliseconds. When series inductance was 
added to the system, the response time slowed to approxi- 
mately five milliseconds. Thus if one decides to add series 
inductance to the motor to reduce txie current ripple 
effects, that decision must be tempered by the fact that the 



transient response cf the system will change due to the 
change in the motor’s electrical time constant. 



B. EEDUCTION OF VELOCITY BIPPLE 



Optimization of motor performance will reguire that the 
ripple content of the motor velocity at steady state be 
minimized. To reduce the ripple, it will then be necessary 
to reduce the speed tolerance settings which will establish 
the pulsewidth modulated input signal to the motor. To 
determine the effects, if any, on motor performance (other 
than the reduction of velocity ripple), it was necessary to 
perform a number of simulations with varied speed tolerance 
settings. 

Three simulation trials were conducted, with the speed 
tolerance set at ±5 rpm, ±1 rpm and ±. 1 rpm. Table IV 
summarizes the results of the simulations. The data for 
Table IV is plotted in Figures 4.4 and 4.5. 



TABLE I? 

Performance Trials for Various Speed Tolerance Settings 



Speed Tolerance (rpm) 

±5 ±1 ±.1 



Load Torque 
( 02 “in) 


“(If'" 


form 

factor 




form 

factor 




f orm 
factor 


05. 0 


2. 19 


2.860 


0.89 


1.765 


.214 


1.7 09 


16. 0 


3.00 


1.819 


1 .16 


1. 138 


.242 


1.121 


32. 0 


3.44 


1.372 


1.22 


1. 033 


.235 


1 .026 


48. 0 


3.57 


1.283 


1.21 


1. 012 


.228 


1.011 


64. 0 


3.71 


1.095 


1.26 


1.005 


.228 


1 .005 


80. 0 


3.79 


1.054 


1.29 


1.003 


.250 


1.002 


It is 


apparent 


that reducing 


the speed 


tolerance does 


indeed reduce the 


ripple 


content 


of the 


motor’s 


speed . 



Additionally, as the tolerance is reduced, so does the form 



factor of the motor(for specific loading) , indicating higher 
motor performance efficiency. Thus as speculated before. 



with lower speed tolerance 
performance. Of course t 
obvious even prior to cond 
to the fact that the gre 
greater the inertia obtain 
rotor before a speed lim 
amounts of energy would h 
increase the motor’s speed 
speed. 

C. COMaZNTS 

I^e have seen where t 
modelled brushless dc moto 
addition of series induct 
speed tolerance settings 
There are reasons to assum 
ments noted in this chapter 
in an operational cruise 
chapter of this report disc 
when reviewing this work an 
in the use of brushless dc 



s, tnere exists improved motor 
his should have been intuitively 
ucting the simulation trials due 
ater the speed tolerances, the 



ed 


by 


the 


rotating mass of the 


it 


is 


r ea 


ched and thus greater 


a ve to 


be 


expended in order to 




back 


to 


the current commanded 



he simu 
r has b 
ance and 
of the 
e that 
may no 
missile 
usses wh 
d prior 
motors i 



lated performance of a 
een improved through the 
the optimizing of the 
chopper control logic, 
the performance enhance- 
t necessarily be realized 
scenario. The final 
ere caution need be taken 
to applying these results 
n aerospace applications. 
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VELOCITY RIPPLE STUDIES 



o 




Figure 4.4 Motor Velocity Ripple vs. Load Torque 
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FORM FACTOR 
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FORM FACTOR STUDIES 




Figure 4.5 Form Factor vs. Load Torque (VTOL varied) 
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V. RECOMHEMDATIONS FOR FORTHER STUDIES 



As stated in Chapter 3, 
missile’s control logic wou 
the motor controller unit, 
for all missile applicati 
generate torque commands (o 
as current is directly prop 
motor to provide a constan 
missile maneuvers. Ihe ef 
simple control scheme studi 
be fairly significant, and 
detail. 

A. BEYOND SPEED CONTROL 

One major change that 
controller if torque comm 
missile control logic is 
loops would necessarily ha 
amplifier. Of course, the 
commands is that the oper 
necessity would have to be 
braking and the regenerat 
The importance inherent wi 
of these modes is best just 
ative braking mode serves 
torque applied to the r 
missile’s supply batteries, 
have logic blocks which cou 
mode would be appropriate 
commands to the motcr cont 



the assumption was made that the 
Id generate velocity commands to 
This may not necessarily be true 
ons. Certain controllers may 
r equivalently, current commands, 
ortional to motor torque) to the 
t torque to affect the required 
feet that this would have on the 
ed in this report could prove to 
will now be looked at in closer 



would 
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requ 


ired of 
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to be 
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from 
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that 
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curr 
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ve to 


be 
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PW!1 
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^ue 
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env 
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motor 


by 


defin 


ed t 
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lude the 


plugging. 



ive braking modes of operation, 
thin the inclusion and modelling 
ified by noting that the regener- 
not just to control the braking 



otor 


but also 


1 to 


rech 


arge the 




The model 


would 


1 then 


need to 


Id 


recognize 


when 


each 


operating 


and 


then generate 


the 


requisite 


rol 


logic . 


The 


ideal 


ized step 
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response of the fin actuato 
position is shown in Figu 
operating modes shown as a 

Since 1) a position loo 
accounting has been made f 
than motoring in the forwa 
recommended that these are 
rately model the motor as 
missile fin actuation syst 
modelling the position cont 
represent the next logical 
dynamics of the electromech 

B. SIMUIATION DEFICIENCIES 

Studies of the modelled 
the response of the syste 
constant load conditions, 
conditions, the system mos 
a series of varied load con 
course towards the target, 
be a function of the aerod 
subjected during flight, 
acceleration, etc. In orde 
under the dynamics of fligh 
tion loop should be follow 
so that the motor’s behavio 
more closely related to its 
Other studies that would 
certainly are not limited 
switching affects on the 
designing the hardware r 
controller logic, as we 
Studying of the effects of 



r system to command chance in fin 
re 5. 1 with the required motor 
function of flap position, 
p has yet to be closed and 2) no 
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oper 


ation 


other 


rd 


direc ti 


on. 
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Of course, 
t undoubta 
ditions as 
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r then to 
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prove wo 
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closing a 



ve been concentrated on 
input commands under 
under normal operating 
bly would be subject to 
the missile steers its 
on the motor would then 
o which the missile is 
sile speed, attitude, 
model the system better 
fort to close the posi- 
precise load studies, 
tudied within a context 
operating environment, 
rthwhile include (but 
e study of the voltage 
r's power transistors, 
implement the motor 
associated software, 
phase-locked servo loop 



Ffn Angle 




Figure 5. 1 Theoretical Fin Actuator Response 
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for velocity control might 
system’s precise speed regu 

One caution worth noti 
hasis for all simulation st 
flux back emf waveform. Xh 
sinusoidal in nature, and 
addition of fundamental and 
the motor parameters which 
velocity. It is recommende 
advanced motor model which 
utilized. 

C. SUMMARY OF RESOIIS 

Pulsewidth modulation 
method of accurately and r 
of a brushless dc motor. F 
power losses in the switchi 
thus allowing for smaller p 
sinking, which translates 
design and implementation, 
tance was shown to a have 
current ripple, reducing i 
simulations conducted witho 
the limit cycle pulsewidth 
a superior method for imp 
Areas where further researc 
presented, including speci 
studies to this thesis. 



also be considered due to the 
lation capaoility. 
ng is that the model used as a 
udies assumed a linear, average 
e back emf signal is in actuality 
would therefore result in the 
harmonic frequency components to 
were studied, such as current and 
d that for further studies a more 
simulates sinusoidal back emf be 



has been shown to be a viable 
eliably controlling the velocity 
orm factor studies indicated that 
ng transistors may be minimized, 
ower transistors and reduced heat 
into reduced costs in controller 
The addition of series induc- 
a very definite impact on motor 
t significantly in comparison to 
ut the series inductance. Use of 
modulation scheme was shown to be 
lamenting pulsewidth modulation, 
h efforts may continue have been 
fic recommendations for follow-on 
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APPENDIX A 
THE MOTOR MODEL 



The motor that was mode 
brushless dc motor. The 
motor have been shown in Fi 
a three phase, four pole m 
accomplished electronically 
effect position sensor dev 
transistors. The switchin 
mechanical rotation. 

The back EMF signal wa 
tional to the motor speed, 
the waveshape was not taken 
What follows is a brief d 
that were added to the ba 
description of program v 
program as well as cbserva 
tion of the simulation prog 



ercially available 
ed curves for the 
.4. The motor is 
commutation being 
set of three Hall 
a set of six switching 
ace every 30 degrees of 

s assumed to be directly propor- 
The actual, sinusoidal nature of 
into account in the model used, 
escription of certain procedures 
sic motor simulation program, a 
ariables added to the initial 
tions made concerning the execu- 
ram . 



lied was a coram 
current and spe 
gures 3.3 and 3 
achine with the 
utilizing a 
ices and 
g takes pi 



A. PEOCEDOEES ADDED TO THE BASIC PROGRAM 

1 • Pro ced ur e I CLIP 

Procedure ICLIP was added to the model to account 
for the addition of the freewheeling diode FWD. Currents 
were then necessarily clipped (negatively going) at zero 
amperes, thus preventing the circulation of negative 
currents. 

2 • Proced ure VC IP 

This procedure takes as an input the velocity error 
signal and yields the motor's input voltage based upon 
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coci manded speed and the est 



within 


this procedure 


that 


system 


is established. 




3. 


Procedure RESET 






As the output 


of t 


rotor 


position counts 


upw 


procedu 


re was required 


to 


degrees 


when 360 degrees 



achieved. EESET keeps trac 
rotor turns past the 360 de 
tion to update the variable 
zero to 360 degrees. The 
ultimately be realized when 
the system. 



B. MOTOB EAEAMETEES 

The parameters which fo 
program to achieve the regu 

1 . VCM D 

VCMD is the command 

2 . ITOL 

VTOL is the velocit 



3. ICLIP 

ICLIP is the motor 
to prevent negative current 

4. THE ST 

THEST is the rotor 
from zero to 360 degrees. 



ablished speed tolerance. It is 
the limit cycle behavior of the 



he integrator block which yields 
ard continuously from zero, a 
reset the rotor position to zero 
of rotation of the rotor was 
k of the number of times that the 
gree point and uses this informa- 
THEST, which cycles ranges from 
importance of this procedure will 
a position loop is closed within 



Ilow are those added to the basic 
ired speed control effects. 

ed motor velocity. 



y tolerance limit in rpm. 



current which has been adjusted 
flow. 



position in degrees which ranges 
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NOTES ON PROGRAM EXECOTION 



A fixed interval integr 
of the variable step Run 
default integration technig 
technique was used as it d 
accurately as the varuable 
50% less computer time. Th 
as . 000001 seconds. 

Since the switching fr 
order of 2500 KHZ, to ac 
during the pulse on and off 
than 50 microseconds (typic 
was required. Because of 
the microscopic detail of 
variances present in the c 
greater than one second, a 
approximately 5500 lines o 
for this report the limita 
major problem, but could 
studies. Of course, as the 
the requirements for such d 
ation may not be present, 
used, thus allowing studies 



ation tech 
ge- Kutta 
ue in CSHP 
emonstra te 
step metho 
e integrat 



nigue was used in lieu 
method supplied as a 
-III. The trape 2 oidal 
d itself to perform as 
ds but used better than 
ion interval was chosen 



equency of the motor was on the 
curately observe motor behavior 
periods a print interval of less 
ally 20 microseconds was chosen) 
this, it was difficult to observe 
motor operation in terms of the 
urrent and speed for periods of 
s the CSMP program is limited to 
f output. For the studies made 
tion encountered did not pose a 
prevent an obstacle to further 
studies of this system advances, 
etailed assessment of motor oper- 
and larger print intervals may be 
to be of greater duration. 
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APPENDIX B 



CSMP SIMULATION PROGRAM 



//ASKCSMP2 JOB ( 1 102, 0 1 1 6J , * AS K IN AS AAA A * , C L ASS = C 

//*MAI^ 0RG=NPGVM1. 1102P 

// EXEC CS^FXV 

//X.COPPRINT DC CUMPy 

//X.SYSPRIM DC CUM>y 

//X.SYSIK DO 

INJTI AL 

CGNSTAM KT = 15.9. BM = 0.00015, BL = 0.0, JL = O.C, . 
N = 1.0, JN = C.OOl, KB = 0.112, TLP = 104. 

PARAMETER LA = .0016, RA = 2.740, ... 

TLF = 05., VTCL = 1., VCMC = 1400. 

* A1 = LA/RA -- THE INVERSE ELECTRICAL TIME CCNSTANT 

* A2 = J/B — THE INVERSE MECHANICAL TIME CONSTANT 

* VTOL — THE VELOCITY LIMIT SETTING 

* VCMO — THE CCMMANDEO MOTCR VELOCITY 

* TLF — LOAD TCBCUE IN OZ-INS 

* KB -- EACK EHf CC^STANT 

* KT — TORCUE CCNSTANT 

* N — GEAR RATIO 

* JM — ROTCR MCMENT OF INERTIA 

=*= TLP — LOAD TCRCUE WHERE PEAK POWER CUT CCCURS 



NO SORT 

BLP = BL/(N*=«2J 
JLP = JL/<N**2I 
J = JM + JLP 
B = EM + ELP 
A1 = 2.0 * LA / RA 
A2 = J / E 
THRST = 0.0 
JFAC =0.0 

DYNAMIC 

VERR = VCMO - WMRFM 
VIF= VCLP * STEPtC.OJ 
TLl = RAMP(O.CJ 
TL2 = RAMP(.Cl) 

TL = TLF * ( Til - TL2J/. 01 

VIB = 0.0 

VI N = VI F + VIE 

VINl = V IN - V£MF 

VIN2 = VIM ♦ (l.O/RAJ 

IM = REALPLI 0.0,A1,VIN2) 

TM = ICLIP * KT 

* TNI = TM - TL 

TN2 = TNI ♦ ( 1.0/EJ 
WM = REALPL iO.O,A2,T N2i 
WMRPH = WH <30. /3. 14159) 

WMRPMR = WMRPH/N 
VEMF = WH * K£ 

THETA = INTGRL<0.C,WMi 

THDEG = THETA ♦ ( 180.0/3.14159265) 



THCON = 7HRS7 
PWR = * TM 



;js:{t j{ :ji j{t ;{!;{; :{j :{5 ;j> :<x;45 

* PROCECURE ICLIP SI^^ULATES PRESCENCE OF * 

* FREEWHEELING LJCDE 6Y CLIPPING THE NEGATIVE « 

* GOING CURRENT AT C AMPS * 

PROCEDLRE ICLIP = CLPCRilMJ 

If UM .LE. O.CJ ICLI P = 0.0 
If UM .GT. 0.0) ICLIP = IM 
ENDPROCEDURE 



* ♦ * jj: 4: * -4 * :{t 4s ^{i * :{t 4 ?!: :{t :{! :<t ;{£ :{t :<t :{t X: * ^it * * 4 * * =.< * ^ ♦ 4> =!t 

* PROCEDURE VCLF ESTABLISHES LIMIT CYCLE BEHAVICR. ^ 

* IF THE VELOCITY IS BELOW THE ESTABLISHED SETPOINT,* 

* THE INPUT VOLTAGE IS PULSED ON ( VCLP = 30) AND IF * 

* THE VELOCITY IS AEJV E VTCL + VCMC, THE INPUT * 

* VOLTAGE JS PULSED OFF {VCLP = 0.0) * 

44 5{ ^ ^ sf£« 4 ********** #* 44 



PROCEDLRE VCLP = VVV(VERR, VTCU 
MVTCL = -l.C * V70L 
IF (VERR .GT. VTCL) VCLP = 20.0 
IF (VERR .LT. MVTOL) VCLP = 0.0 
ENDPROCEDURE 

444444 4444444 4*444:4444 4444 444 4 4 4444444 44444 4444 444 

* TNI SETS THE SIGN OF THE LOAD TORQUE TO ENSURE * 

* THAT THE TORCLE OPPOSES THE MOTOR TORQUE WHEN * 

* THE NCTOF OPERATES IN THE FORWARD DIRECTION * 

* AND ADDS TO THE MOTOR TORQUE WHEN * 

* THE MOTOR OPERATES IN THE REVERSE DIRECTION * 

44 444 4 4444 4444 4 444444 4 4444 4 44 4 444 4 444 4 4 44 4 4 444 4 444 



PROCEDURE TNl = f WCEWC( V IN,TM ,TL) 

IF( V IN.L T.O.C ) GO TO 10 
TNI = TM - TL 
GO TO 15 

10 TNI = TM + TL 
15 CONTINUE 
ENDPROCEDURE 

44 4444 4444 444 4 4444444 4 4444 444 4 444 4444 4 444 4 4 444 4 444 
4 PWRA HAS BEEN ADDED TO THE PROGRAM TO ACCCUNT 4 
4 FOR THE NONLINEAR ASPECT OF THE MODELLED 4 

4 MOTOR'S OUTPUT POWER FUNCTION 4 

:{.«:«£# jjt:>4s4 4 4 444 44 4 4444 4 44 4 4444444 4 44444 444 4 444 



PROCEDURE PWRA=T«IDOL( PWR, TL, TLP) 
IF(TL.GE.TLP) GO TO 20 
PWRA = PWR 
GO TO 25 

20 TwID = 4.50 4 (TL - TLP) 
PWREX = EXP(-TWID) 

PWRA = PWR 4 PWREX 
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25 CONTINUE 
ENDPROCEDURE 



5{:3{x ;{*# :Ji :f:{4 :{t :{t ;{t 

« THIS PROCEDURE h.AS ADDED TO RESET THE ROTOR'S 

* POSITION AFTEf IT REACHES 360 DEGREES 

* BACK TO 0 DEGREES 

PROCEDURE THRST=RES£T( JFAC ,THDEG) 

TS = JFAC ♦ 360.0 
THRST = THOE6- TS 
IF<THRST.LT.26C.0i GO TO 40 
JF/SC = JFAC ■* 1 .0 
40 CONTINUE 
ENDPRGCEDURE 



TERMINAL 

TITLE BASIC CC MOTOR 
TIMER FINTIM = .015, 
PRDEL = .000C4, CELT 
METHCD TRAP2 
* OUTPUT WMRPM, THETA, 
PRINT ICLIP,^CLP, 
LABEL MOTCR SPEED 
PAGE MERGE 



SYSTEM 
OUT CEL = , 
= .000001 



WMRPM, 
DOE TO 



THRST 

THRST 
STEF J 



PAGE XYPLOT 
END 
STOP 



ENDJOB 



00004, 

NPUT 



• • • 
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APPENDIX C 



DATA ANALYSIS PROGRAM 



SENTRV 

$JOB CLRRtXREF 

C 

C THIS PROGRAM I' DESIGNED TO YIELD THE AVERAGE CURRENT 
C AND THE RMS CURRENT AS WELL AS K-FACTOR FCR GIVEN 
C SIMUIATICN OF / CCMPUTER MODEL OF A PwM 
C CONTROLLED BRUSHLESS DC MOTOR 
C 

REAL 10, ICC, RES, I NCU C ,RPM AV E , PI , T IME , T IM E 2, FR EU 
REAL EO , IAV,12RMS,A,B,ZN,KFAC,IRMS,iCN,ICFF,IM,IN 
REAL VIN,TC,P,TL,DC,PW,iC£N 
FI = 3.1^159 

RES IS THE RESISTANCE OF THE COILS 
INOUC IS THE INDUCTANCE BETWEEN THE TWO TERMINALS 
VIN IS THE INPUT VOLTAGE CF THE MOTOR 

RES = 3.2^C 
INDJC = .C013 
VIN = 30. C 

FREC = FREQUENCY. OF PwM wAV E SHAPES 
DC REFERS TC DOlV CYCLE CR ‘ON TIME* OF THE WAVE 

TL IS THE LCAD TCRQUE FCR THE GIVEN RUN 
10 IS THE PEAK CURRENT ACHIEVED 
100 IS THE MINIFIUM CURRENT SEEN BY THE MOTOR 
TO IS THE AMCLNT OF TIME THAT THE CURRENT WAVESHAPE 
IS LESS THAN ZERO 

FREQ = 1U3. 

DC = .65^ 

TL = 9.6 
RPMAVE - 1400 
10 = 4.47 
ICC = 3.42 
TO = .OOC 

ZN IS THE SYSTEM ELECTRICAL TIME CONSTANT 

ZN = INDLC/RES 
PW = (l.O/FREC) * DC/100. 

TIME = PW 

TIME2 = 1. C/FREQ - PW 

EO IS THE BACK EMF 

EC = 1.631 * RPMAVE *.0558 *PI/3C.O 

IN IS THE THECRETICAL MAXIMUM OF THE CURRENT 
IM IS THE THECRETICAL MI MMUM OF THE CURRENT 
A IS THE PULSE CN TIME 
B IS THE PULSE CFF TIME 
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IW = £0/P£S 
JN = (VI^ - ECJ/RES 
C 

A = Zh * /SLCGUIN - IGOJ/dN - ICd 
E = 2N * ;iLCGi(IO + I^i/dOO + IMd 



F = A + E TC 

lAV = ( (A<IM - (E*IW JJ / P 

IDEN = Z^*(IC - IOOJ<'(JK + IMJJ/P 

I2RMS = ((A*IN*IN + - IDEM 

IRMS = SCFT(I2RMSJ 

KfAC = IFNE/IAV 



i«RIT E ( 6 ,lCi 

10 FCRMATdA .'STATISTICS FOR PwN CONTROLLED CC MCTORi'J 
VRITE (6,iO)FPEQ 

20 FCRMAT( !>, 'FREQUENCE = ' ,F8 .2 H Z ' , ) 

FRITE (6,22iTL 

22 FORMAT! 1> t' LOAD TORQCE = 'fF8.2i 
i^RITE (6,50) IN, IM 

50 FORMAT (1>,'JN: ',F10.3,' IM: ',F10.3) 

WRITE ( 6, ICC) lA VfIRMS 

100 FCPMATd X lAV = ',F8.4,/,' IRMS = ',F8.4J 
WRITE ( 6,iC03 KF AC 

200 FORMAT ( IX , ' K FACTOR ^ ',F8.5) 

STOP 

END 

SENTRY 
$E NTRY 
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APPENDIX D 



SAMPLE OUTPUT FOR CSMP SIMULATION 



TIM E 


ICLIP 


VCLP 




.0 


.0 


30.000 


.0 


4.0C000D-05 


.36860 


30.000 


1.1220 


8.0C0C0C-03 


. 7 244 6 


30.000 


4.4361 


1.2C000C-04 


1.0677 


30.000 


9.8656 


l,6C000i:-04 


1.3984 


30.000 


17. 334 


2 .0 0000 C-0 4 


1 .7168 


30.000 


26. 765 


2.4COOOD-04 


2.02 3 0 


30.000 


38. 086 


2.8C000C-04 


2.31 71 


30.000 


51. 221 


3.20000C-04 


2.599 3 


30.000 


66.099 


3.6C000D-04 


2. 669 7 


30.000 


82.647 


4.0C000C-04 


3. 1285 


30.000 


100. 79 


4.4C000C-04 


3.3759 


30.000 


120 .47 


4.8C000D-04 


3.6119 


30.000 


141 .61 


5 .20000 C-0 4 


3.836 9 


30.000 


164. 13 


5.6CCOOD-04 


4.0510 


30.000 


187 .99 


6.0C000 C-0 4 


4 .2544 


30.000 


213 . 10 


6.4C000C-04 


4.447 1 


30.000 


239 . 40 


6. 8C000D-04 


4 . 62 9 5 


30.000 


266 . 84 


7.2C000C-04 


4.6018 


30.000 


295.35 


7.6COOOC-04 


4.9641 


30.000 


324.86 


8.0C0C0C-04 


5.116 5 


30.000 


355.32 


8.40000C-04 


5.2594 


30.000 


3 8o . 68 


8.8C000D-04 


5.3929 


30.000 


418 .86 


9.20000C-04 


3.5172 


30.000 


451.82 


9.60000C-04 


5 .632 5 


30.000 


485 .49 


l.OCOOOD-03 


5.739 1 


30.000 


519 .84 


1 .O-^OOOC-03 


5.6370 


30.000 


554. 79 


1.0 6000 C-03 


5.926 6 


30.000 


590 .31 


1. liOOOD-03 


6.C080 


30.000 


626.34 


1.16000C-03 


6.0814 


30.000 


662.83 


1.2C0OOO-O3 


6. 147 1 


30.000 


699 .74 


1.2^0000-03 


6.2052 


30.000 


737. Cl 


1.2 6000C-03 


6.2560 


30.000 


774.61 


1. 320000-03 


6.2996 


30.000 


812.48 


1 .3 6000 C-0 3 


6.33o2 


30.000 


850.59 


1.4COOOD-03 


6.3662 


30.000 


888 .89 


1.4^0000-03 


6.3895 


30.000 


927.35 


1 .4 6000 0-0 3 


6 .406 6 


30.000 


965 .93 


1.520000-03 


6.4175 


30.000 


100 4 .6 


1.560C0C-03 


6.422 5 


30.000 


1043.3 


1.6C0C0C-03 


6.421 7 


30.000 


108 2.0 


1.6^0000-03 


6.4154 


30.000 


1120.6 


1 .6 8000 C-0 3 


6.4037 


30.000 


1159.2 


1.720000-03 


6.3868 


30.000 


1197.7 


1.76000C-03 


6.3650 


30.000 


123 6.1 


1.8 0000 C-0 3 


6.3384 


30.000 


1274.4 


1.8^0000-03 


t.3072 


30.000 


1312.4 


1.86000C-03 


6.2716 


30.000 


1350.2 


1. 92000C-03 


6.231 7 


30.000 


1387 .8 


1.960000-03 


5.9513 


.0 


1424.7 


2 .00000 C-03 


5.5431 


.0 


1459.2 


2.040000-03 


5.143 8 


.0 


1491.3 


2.08000C-03 


4.7534 


.0 


1520.9 


2. 12000 C-0 3 


4.3722 


.0 


1548 .2 
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2 . l ^ OOOC -03 


4 . C 000 


.0 


1572.2 


2 . 2 C 000 C -03 


3.6369 


.0 


1596.0 


2 . 2 ^ 0000-03 


3.2829 


.0 


1016.5 


2 . 260 C 0 C -03 


2.9380 


.0 


1635.0 


2 . 32000 C -03 


2.602 2 


.0 


1651.3 


2 . 3 dOOOC -03 


2.2755 


.0 


1605.7 


2 . 400000-03 


1.9579 


.0 


1678.1 


2 . 4 ^ 0000-03 


1.6493 


.0 


168 8 .5 


2 . 460000-03 


1.2497 


.0 


1697.2 


2 . 520000-03 


1.0591 


.0 


1704 .0 


2 . 560000-03 


.77751 


.0 


1709.1 


2 . 600000-03 


.50474 


.0 


1712.4 


2 . 6 ^ 0000-03 


.2408 0 


.0 


1714.2 


2 . 6 € 0000-03 


.0 


.0 


1714.3 


2.7 2000 0-0 3 


.0 


.0 


1713.8 


2 . 760000-03 


.0 


.0 


1713 .3 


2 . 800000-03 


..0 


.0 


1712.7 


2 . 840000-03 


.0 


.0 


1712.2 


2 . 880000-03 


.0 


.0 


171 1.6 


2 . 920000-03 


.0 


.0 


1711.1 


2 . S 60 COC -03 


.0 


.0 


1710.5 


3 . 0 C 00 O 0 - O 3 


.0 


.0 


1709.9 


3 . 040000-03 


.0 


.0 


1709.3 


3 . O 80 COO - O 3 


.0 


.0 


1708.7 


3 . 120000-03 


.0 


.0 


1708.1 


3.16000 0-03 


.0 


.0 


170 7.5 


3 . 200000-03 


.0 


.0 


1706.9 


3 . 240000-03 


..0 


.0 


1706.3 


3 . 28 OCO 0 - O 3 


.0 


.0 


1705.6 


3 . 320 C 0 C -03 


.0 


.0 


1705.0 


3.36000 0-03 


.0 


.0 


1704.3 


3 . 4 C 0 C 00-03 


.0 


.0 


1703.7 


3 . 440000-03 


• c 


.0 


1703.0 


3 . 480000-03 


.0 


.0 


170 2.3 


3 . 520000-03 


.0 


.0 


170 1.0 


3 . 560000-03 


.0 


.0 


1701.0 


3 . 600000-03 


.0 


.0 


1700.3 


3 . 640000-03 


.0 


.0 


1699.6 


3.68000 0-03 


.0 


.0 


1698.8 


3 . 720000-03 


.0 


.0 


1698.1 


3 . 760000-03 


.0 


.0 


1697.4 


3.8 0000 0-0 3 


.0 


.0 


1696 .7 


3 . 840000-03 


.0 


.0 


1695.9 


3 . 880000-03 


.0 


.0 


1695.2 


3 . 920000-03 


.0 


.0 


1694.4 


3 . 960000-03 


.0 


.0 


1693.0 


4 .00000 0-03 


.0 


.0 


1692.9 


4 . 040000-03 


.0 


.0 


1692.1 


4 . 080000-03 


..0 


.0 


1691.3 


4 . 1 2000 0-03 


.0 


.0 


1690 .5 


4 . 160000-03 


.0 


.0 


168 9.7 


4 . 200000-03 


.0 


.0 


1688.9 


4 . 240000-03 


.0 


.0 


1688.1 


4 . 280000-03 


.0 


.0 


1687.3 


4 . 320000-03 


.0 


.0 


1686.4 


4 . 360000-03 


.0 


.0 


1685 .6 


4 . 400000-03 


.0 


.0 


1684.7 


4.4 4 0 000-03 


.0 


.0 


1682.9 


4 . 480000-03 


.0 


.0 


1683.0 


4 . 520000-03 


.0 


.0 


1682.1 


4 . 560000-03 


.0 


.0 


168 1.3 


4 . 6 C 0 COO - O 3 


.0 


.0 


168 C .4 


4.64000 0-03 


.0 


.0 


1679 .5 


4 . 680000-03 


.0 


.0 


lo 78.6 


4 . 720000-03 


.0 


.0 


1677.7 


4.7 6000 0-0 3 


.0 


.0 


1676.7 


4 . 8 C 0000-03 


.0 


.0 


1675.8 
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4 . 8 ^ 000 C -03 


.0 


.0 


1674 .9 


4 , 880 COC -03 


.0 


.0 


1673.9 


4,92000 C -03 


.0 


.0 


1673.0 


4 . 960 C 0 D -03 


.0 


.0 


167 2.0 


5 . 0 C 000 C - O 3 


.0 


.0 


1671.1 


5 . 04000 C -03 


.0 


.0 


1670. 1 


3 . 080000-03 


.0 


.0 


1669 .1 


5 . 120 C 0 C -03 


..0 


.0 


1668.1 


5 . 180 COO - O 3 


.0 


.0 


1667 .1 


5 . 2 COOOC -03 


.0 


.0 


1668.1 


5 .2 4000 0-0 3 


.0 


.0 


1665.1 


5 . 280 COD -03 


.0 


.0 


166 4.1 


5 . 320000-03 


.0 


.0 


1663.1 


5 . 380000-03 


.0 


.0 


1662.1 


5 . 4 COOOO -03 


.0 


.0 


166 1 .0 


5 . 4 ^ 0000-03 


..0 


.0 


1660 .0 


5 . 480000-03 


.0 


.0 


1658 .9 


5 . 520 C 0 C -03 


.0 


.0 


1657.9 


5 . 560000-03 


.0 


.0 


1656.8 


5 . 6 CCC 00-03 


.0 


.0 


165 5.7 


5 . 6 ^ 0000-03 


.0 


.0 


1654.6 


5 . 680000-03 


.0 


.0 


1653.5 


5 . 720000-03 


.0 


.0 


1652.4 


5 . 780000-03 


..0 


.0 


1651.3 


5 . 8 COOOO -03 


.0 


.0 


1650.2 


5 . 8 ^ 0000-03 


..0 


.0 


1649.1 


5 . 880000-03 


.0 


.0 


1647.9 


5 . 920000-03 


.0 


.0 


1646 .8 


5 . 98 CCOO -03 


.0 


.0 


1645.7 


6 . 0 COCOO -03 


.0 


.0 


1644.5 


6 . 0 ^ 0000-03 


.0 


.0 


1643.3 


6.08000 0-03 


.0 


.0 


1642.2 


6 . 120 C 0 D -03 


.0 


.0 


1641.0 


6 . 180000-03 


.0 


.0 


163 9.8 


6 . 200000-03 


.0 


.0 


1638.6 


6 . 2 ^ 0000-03 


.0 


.0 


163 7 .4 


6 . 280000-03 


..0 


.0 


1638.2 


6 . 320000-03 


.0 


.0 


163 5 .0 


6 . 380000-03 


.0 


.0 


1633.8 


6.40000 0-03 


.0 


.0 


1632.5 


6 . 440000-03 


.0 


.0 


163 1.3 


6 . 48000 0-03 


..0 


.0 


1630.1 


6.5 2000 0-03 


.0 


.0 


1628.8 


6 . 580000-03 


.0 


.0 


162 7.5 


6 . 6 C 0 C 0 C -03 


.0 


.0 


1628.3 


6 . 640000-03 


.0 


.0 


162 5.0 


6 . 68 CC 00-03 


.0 


.0 


lc 23.7 


6.7 2000 0-0 3 


.0 


.0 


1622.4 


6 . 78 CC 00-03 


.0 


.0 


1621.1 


6 . 800000-03 


.0 


.0 


161 9.8 


6 .84000 0-03 


.0 


.0 


1618.5 


6 . 880000-03 


.0 


.0 


1617.2 


6 . 920000-03 


.0 


.0 


1615.8 


6 . 980000-03 


.0 


.0 


1614.5 


7 . 0 C 0000-03 


.0 


.0 


1613.2 


7.04000 0-03 


.0 


.0 


1611.8 


7 . 080000-03 


.0 


.0 


1610.4 


7 . 120000-03 


.0 


.0 


1609.1 


7 . 160000-03 


.0 


.0 


160 7 .7 


7 . 2 COOOO -03 


.0 


.0 


160 6.3 


7 . 240000-03 


.0 


.0 


1604.9 


7 . 280000-03 


.0 


.0 


160 3.5 


7 . 320000-03 


.0 


.0 


160 2.1 


7 . 360000-03 


.0 


.0 


1600.7 


7 . 4 COOOO -03 


.0 


.0 


1599.3 


7 . 440000-03 


.0 


.0 


1597.9 


7.4 8000 0-0 3 


.0 


.0 


1598.4 



60 



7 . 52000 C -03 


.0 


.0 


1595 .0 


7 . 5 ^ 000 C -03 


.0 


.0 


1593.5 


7 . 6 C 0000-03 


.0 


.0 


1592.1 


7 . 64 CC 0 C -03 


.0 


.0 


1590 . o 


7 . 6 iC 00 C -03 


.0 


.0 


1569.1 


7 . 72 C 00 C -03 


.0 


.0 


1587.6 


7 . 76000 C -03 


.0 


.0 


1586.1 


7 . 8 C 0000-03 


.0 


.0 


1384 .6 


7 . 8 ^ 000 C -03 


.0 


.0 


1583.1 


7 . 88000 C -03 


.0 


.0 


158 1 .6 


7 . 92 C 00 C -03 


.0 


.0 


1580.1 


7.96000 C -03 


.0 


.0 


1578.6 


8 . 0 C 0 C 0 D -03 


.0 


.0 


157 7.0 


8 . 0 ^ 000 C -03 


.0 


.0 


1575.5 


8 . 08000 C -03 


.0 


.0 


1574.0 


8 . 120000-03 


.0 


.0 


1572.4 


8 . 160 C 0 C -03 


.0 


.0 


1570.8 


8 . 2 C 000 C -03 


.0 


.0 


1569 .3 


8 . 2 ^ 0000-03 


.0 


.0 


1567 .7 


8 . 280000-03 


.0 


.0 


1566.1 


8 . 320000-03 


.0 


.0 




L 564.5 


8 . 360000-03 


.0 


.0 




L 562.9 


8.40000 0-03 


.0 


.0 




L 56 1.3 


8 . 440000-03 


.0 


.0 




L 55 S .6 


8 . 460000-03 


.0 


.0 




L 558.0 


8 . 520000-03 


.0 


.0 




L 556 .4 


8 . 560000-03 


.0 


.0 




L 554.7 


8.6 0000 0-0 3 


.0 


.0 




1553 .1 


8 . 640000-03 


.0 


.0 




L 55 1.4 


8 . 660000-03 


..0 


.0 




L 549.8 


8 . 720000-03 


.0 


.0 


J 


L 548.1 


8 . 760000-03 


.0 


.0 




546 .4 


8 . 800000-03 


.0 


.0 




L 544.7 


8 . 840000-03 


.0 


.0 




L 543.0 


8 . 860000-03 


.0 


.0 




L 541 .3 


8.9 2000 0-03 


..0 


.0 




1539 .6 


8 . 960000-03 


.0 


.0 




L 53 7.9 


9 . 0 C 000 C -03 


.0 


.0 




L 53 6.2 


9.04000 0-0 3 


.0 


.0 




.53 4 .4 


9 . 060000-03 


.0 


.0 




L 53 2.7 


9. 12000 0-03 


.0 


.0 




L 53 C .9 


9 . 1 60000-03 


.0 


.0 
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